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Abstract 



Using Ax 2 which is defined by the difference of the number of events with 
the CP phase 5 and the hypothetical one with 5 = 0, we discuss correlations 
of errors of the CP phase and other oscillation parameters as well as the 
matter effect in measurements at neutrino factories. By varying the oscillation 
parameters and the normalization of the matter effect, we evaluated the data 
size required to reject a hypothesis with 5 = at 3aCL. The optimum muon 
energy and the baseline depends on the magnitude of #13, the background 
fraction, the uncertainty of the normalization of the matter effect, but in 
general lie in the ranges 20GeV < < 50GeV, 1000km< L <3000km. If 
we assume that the uncertainty of the matter effect is as large as 20% then 
the optimum values may be modified to E^ < lOGeV, L < 1000km due to 
the strong correlation of 5 and the matter effect. We show analytically that 
sensitivity to CP violation is lost for E^ <C lOGeV or for E^ » 50GeV. We 
also discuss the possibility of measuring CP violation at the upgraded JHF 
experiment by taking all the error correlations into account, and show that it 
is possible to demonstrate 8 7^ at 3<rCL for #13 > 3°. 
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I. INTRODUCTION 



There have been several experiments J1|-|TT| which suggest neutrino oscillations. 
It has been shown in the two flavor framework that the solar neutrino deficit can 
be explained by neutrino oscillation with the set of parameters (Am|, sin 2 29 & ) ~ 
(£>(lO- 5 eV 2 ),0(l(r 2 )) (small angle MSW solution), (£>(10- 5 eV 2 ), 0(1)) (large angle MSW 
solution), or (£>(10- 10 eV 2 ), 0(1)) (vacuum oscillation solution), and the atmospheric neu- 
trino anomaly can be accounted for by (AmL m , sin 2 26* a ^ m ) ~ (10~ 2 ' 5 eV 2 , 1.0). 

In the three flavor framework there are two independent mass squared differences and 
it is usually assumed that these two mass differences correspond to Am^ and Am^ m . 
Throughout this paper we will assume three neutrino species which can account for only the 
solar neutrino deficit and the atmospheric neutrino anomaly [J Without loss of generality 



we assume |Am 21 | < |Am 2 9 | < 



°32\ 



| Am 2 x | 



related to the mass eigenstates by U ( 
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With the mass hierarchy jAm^il jAm^l there are two possible mass patterns which 
are depicted in Fig. la and lb, depending on whether Am 2 2 is positive or negative. If the 
matter effect is relevant then the sign of Am 2 2 can be determined by distinguishing neutrinos 
and anti-neutrinos. The Superkamiokande experiment uses water Cherenkov detectors and 
events of neutrinos and anti-neutrinos are unfortunately indistinguishable, so the sign of 
Am 2 2 is unknown to date. 



It has been shown in the three flavor framework [13-15 



of the CHOOZ reactor data 
| n| and the atmospheric neutrino data implies very small 813, i.e., sin 2 26 13 < 0.1 which is 



essentially the result of the CHOOZ data. When | 1 is small, the MNS matrix looks like 
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where 6*12, 6*23 have been replaced by 9 & and respectively. 

Measurements of #13 and even CP phase S may be possible in the near future projects 
|T7| . [T8|| . Neutrino factories HT9| , ^0| have been proposed to measure 9 13 , the sign of Am| 2 and 
CP phase 5, and a lot of people have worked out many aspects in detail p]|-p3|. 

Since measurement of CP violation is the main motivation for neutrino factories, it is 
important to examine errors of the CP phase S taking into consideration all possible sources 



1 To explain the LSND anomaly [11] one needs at least four neutrino species. 
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of uncertainties in the oscillation parameters as well as matter effects. In this paper we 
discuss correlations of errors of 5 and the matter effect, as well as Qm and Am 2 .^. We 
find that the correlation of 8 and the matter effect is strong for (£^=50 GeV, L ~3000km), 
which has been advocated as the best choice, and if the uncertainty of the matter effect is 
larger than 5 % then baselines shorter than 3000km are advantageous. 



There have been works p2| , p4| , p5| which advocated the advantage of conventional low 
energy neutrino beams over neutrino factories with high energy, and there does not seem 
to be consensus among the community of neutrino physicists about the optimum neutrino 
energy and the baseline. One of the purposes of this paper is to show that neutrino fac- 
tories with high or medium muon energy (20GeV< <50GeV) are more advantageous 
over experiments with low energy (£ M <ClOGeV). We evaluate the statistical significance of 
possible CP violation in neutrino factory type experiments for a range of the muon energy 
0.5GeV< E^ <50GeV and the baseline 10km < L < 10 4 km and we show that the case £ M < 
a few GeV is always inferior to the option with higher energy. On the other hand, Lipari | 3H[ 
has shown that long baseline experiments lose the sensitivity for CP violation for extremely 
high energy. We will show analytically that the two kinds of A% 2 that we will use, one of 
which is defined through the difference of 5 ^ and 5 = and the other of which through 
T violation, decrease for large muon energy E^ ^>50GeV. 

We also analyze the possibility to measure CP violation at the upgraded JHF experiment 
with 4MW power and a one mega ton detector. We show that measurement of CP violation 
is possible for a relatively large value of #13. It is possible mainly because the detector size 
is so huge. 

In section 2 we analyze numerically the effect of the CP phase 5 in neutrino factories 
by introducing A% 2 which is defined as the square of the difference of the number of events 
with 5^0 and that with 5 = 0, divided by errors. We show that the optimal muon 
energy and baseline vary depending on the magnitude of #13, the background fraction and 
the uncertainty of the matter effect. In section 3 we discuss T violation assuming ideal 
polarization. In section 4 we present analytical treatments of the two kinds of A% 2 introduced 
in sections 2 and 3, and show that these A% 2 decrease as the muon energy becomes very 
small or very large. In section 5 we discuss the possibility of measurement of CP violation 
at the upgraded JHF experiment. In section 6 we summarize our results and we briefly 
comment on the uncertainty of the matter effects from the viewpoint of geophysics. 



II. ANALYSIS OF CP VIOLATION 



2 Some of the correlations of errors of 5 and other oscillation parameters have been discussed in 
the past [22, 29 1, but they used relatively optimistic sets of the oscillation parameters, and the 
behaviors we obtain with the most up to date best fit values for the oscillation parameters are 
sometimes different from theirs. 
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A. Definition of Ax 2 



Our strategy here is to examine whether a hypothesis with a vanishing CP phase is 
rejected or not by taking into consideration all channels v e — > z/ M , v e — > u^, —>■ and 

— > *V F° r this purpose we define Ax 2 to test a hypothesis that a CP phase is given by 
5 in the case where the true value is 5: 
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where i?^ is the muon energy, L is the length of the neutrino path, iVo is the number of 
the target nucleons times the number of useful decays of muons, ct u n(E u ) and a P ]y(Ep) are 
the (anti-)neutrino nucleon cross sections. We adopt the cross section which is the sum 
of those |J7| of the quasi elastic scattering, one pion production, and inelastic scattering, 
where double counting of the latter two is suitably subtracted |38 |. Throughout this paper 
the threshold energy is assumed to be O.lGeV which is close to what has been assumed for 
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liquid argon detectors [2{| and which may be realized in possible mega ton water Cherenkov 



detectors 39.40 



The number of the free parameters in the present case is six (5, # i2 , #13, # 23 , Am 



2 

21) 



Am^), but the density N e (x) = Y e (x)p(x) of electrons is not known exactly (Y e (x) is the 
ratio of the number of electrons to that of protons and neutrons, and p(x) is the density 
of the Earth at a distance x from the beam production point), so we have to vary N e (x) 
also. Here for simplicity we assume the PREM (Preliminary Reference Earth Model) 
and vary the overall normalization of the PREM: 



A(x) = CA (x) = y/2CG F Y e (x)N e 



x) 



where C = 1 corresponds to the PREM. We have to consider correlations of errors of the 
CP phase and six other quantities and taking into account all these errors we obtain the 
probability of rejecting a hypothesis 5 = 0. To do that we look for the minimum value 
of Ax 2 (0m, Amy, 5, C; 9 k e, Am^, 5, C) by varying the six parameters (# i2 , #13, # 2 3, Am^, 
H, C): 



A X 2 min = min_ A X 2 (6 ke} Amt e} 5,C;6 ke} Aml l ,5 = 0,C), 
where C stands for the overall normalization of the electron density. 



B. Correlations of errors of 5 and other parameters 

Let us first discuss correlations of two variables (5, X) where a parameter X stands for 
C, #13, #12, #23, Am 21 and Am 2 2 . 

We have studied numerically correlations of errors between 5 and the other os- 
cillation parameters (9 k £, Amy) as well as the normalization C of the matter ef- 
fect for the case where the central values for these parameters are those of the 
best fit point, i.e., sin 2 29 12 = 0.75, Am^ = 3.2 x i(T 5 eV 2 ; sin 2 2# 23 = 
1.0, Am 2 2 = 3.2 x 10~ 3 eV 2 , C=1.0 and we have used a reference value 
8°. The values of A X 2 (9 12 , #13, # 23 , Am 2 1? Am 2 2 , 5, C; 9 12 , #13, #23, Am 2 1; Am 2 2 , 5, C), 



Ax 2 (#i2, #13, #23, Am 2 !, Am 2 2 , 5, C 
Ax 2 (#i 2 , #13, #23, Am 2 1; Am 2 2 , 5, C 
Ax 2 (#i 2 , #13, #23, Am 2 1; Am 2 2 , 5, C 
Ax 2 (#i2, #13, #23, Am 2 1; Am 2 2 , 5, C 
Ax 2 (#i2, #13, #23, Am 2 1; Am| 2 , 5, C 



#12, #13, #23, Am 2 1? Am 2 2 , 5_, C), 
#12, #13, #23, Am 2 x , Am 2 2 , 5_ } C), 
#12, #13, #23, Am 2 1 ; Am§2, 5_ } C), 
#12, #13, #23, Am^, Am§ 2 , 5_, C), 

#12, #13, #23, Am 21 , Am 2 2 , 5, C) are plotted in Fig. 1-6 
in the case of 5 = vr/2 for ^=3, 20, 50 GeV, L=100km, 1000km, 2500km, 6300km, where 
the data size 10 21 /i-10kt is used as a reference value and no backgrounds are assumed. Since 
the number of degrees of freedom is 2, Ax 2 =0.18, 0.34, 0.73 correspond to la, 90%, 99% 
confidence level to reject a hypothesis with 5 = 0. 

As can be seen in Fig. 1, the correlation (5, C) for L ~3000km is strong for #13 = 8°. 
The correlation (5, C) turns out to be small for larger values of Am 21 or for smaller value 
of #13 (i.e., #13 ^5 3°; See Fig. lb), as the gradient of the ellipse in the (5, C) plane 
becomes smaller for larger values of Am 21 . This is why strong correlations were not found 
in p9[ where the set of parameters (sin 2 2# i2 = 1.0, Am 21 = 1.0 x 10~ 4 eV 2 ; sin 2 2# 23 = 1.0, 
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Am 2 32 = 3.5(5,7) x l(T 3 eV 2 ) and sin 2 26 13 = 0.05, = 30GeV were used. If we assume 
that the uncertainty in the overall normalization C is at most 5%, then the correlation (5, 
C) is not so serious, but if we assume that the uncertainty is as large as 20 % then the set 
of the parameters (E^ ~ 50GeV, L ~ 3000km) is not a good option. We will discuss this 
issue later. 

From Figs. 2 - 6, we see that the correlations of (5, O^) and (5, Amy) are not 
large for L > 1000km, E^ >20GeV. As we will show analytically later, the value of 
Ax 2 (#h, Amy, 5, C\ 9m, Am^, 5 = 0,C) increases for E^ > 50GeV and L <C 1000km unless 
we minimize it with respect to 6u and Amy, but because of strong correlations in (<f, # 13 ), 
(5, 623) and (5, Am|), the value of Ax^^, which is minimized with respect to 6te and 
Amy, decreases for ^> 50GeV. 



C. Data size to reject a hypothesis with 5 = 



The quantity Ax 2 - can be regarded as the deviation of x 2 from the best fit point (the 

C for which 
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best fit point in eq. ([!]) is of course 6m = 6ke, Amy 

we have Ay 2 • = 0) and for six degrees of freedom the value of Ay 2 • which corresponds 
mm ' mm 

to 3cr (4cr ) is 20.1 (28.9). From this we can estimate the necessary data size D to reject a 
hypothesis 5 = at 3cr by dividing 20.1 by A% 2 • for each value of 5. On the other hand, 
it is important to include the effect of the backgrounds in the analysis [p^ , p5| ,pE2]] . Here we 
assume that the fraction fs of backgrounds to right sign muon events is given by fs = 10 -3 
or 10~ 5 and that the systematic error of backgrounds is Ob = 0.1 as in for simplicity. 
We also assume the number of muons 10 21 /i- lOkt as a reference value. Thus Ax 2 is modified 
as 
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son statistical |fl2]j . Then we minimize Ax 2 with respect 6^, Arn 2 kl and (7: 
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where the values of oscillation parameters we use in (|3|) are the best fit values in the analyses 
of the solar and atmospheric neutrinos [43] as in Figs. 1-6, and we take #13= 1°, 5°, 8° 
and 5 = tt/2 as a reference value. In varying the overall normalization C we assume 0.95 < 
C < 1.05. We will mention the results for |AC| < 0.1 and for |AC| < 0.2 later. For other 
oscillation parameters, we vary (812, Am^) and (623, Am| 2 ) within the allowed region at 
90%CL of the solar and the atmospheric neutrino data, i.e., 25° < 612 < 41°, 35° < 623 < 55°, 
1.5 x 10~ 5 eV 2 < Am 2 21 < 2.2 x 10~ 4 eV 2 , 1.6 x 10~ 3 eV 2 < Am 2 2 < 4 x 10~ 3 eV 2 . It should 
be emphasized that in minimizing A% 2 in (|3|) all the six parameters are varied at the same 
time, unlike in Figs. 1-6 which are obtained by varying only one of 6ke, Am^, C. 

The result is given in Fig. 7 for a neutrino factory with 0.5GeV < < 50GeV, 10km 
< L < 10000km and for three values of #13= 1°, 5°, 8° and two different values of the 
background fraction fs = 10 -5 , 10~ 3 . The behavior of the figures change a little depending 
on the value of 6*13. For fs = 10 -3 , the sensitivity to CP violation, i.e., the ability to 
reject a hypothesis with 5 = is not optimized by the set of parameters (E^, L) ~ (50GeV, 
3000km), which has been advocated as the best choice, but rather by (E^, L) ~ (20GeV, 
2000km) . This is because with a nonnegligible fraction the contribution of the systematic 

uncertainty ob/bA^- to the total error becomes so large for high energy such as E^ ~ 
50GeV and sensitivity to CP violation is lost. For Jb = 10~ 5 and #13= 1°, on the other 

hand, the contribution of cr B / s ArJ^^ is not so large and the sensitivity is optimized by 
(E^, L) ~ (50GeV, 3000km). We note in passing that we have also optimized the sensitivity 
with respect to the number of energy bins but the conclusion does not depend very much 
on the number of energy bins. This result disagrees with the claim in ]32] . 

We have also evaluated the data size assuming a larger uncertainty of the matter effect, 
i.e., |AC| < 0.1 and |AC| < 0.2. The results for #13 = 8° are shown in Fig. 8. If we have 
to assume an uncertainty of the matter effect which is as large as 20%, then the optimum 
baseline and muon energy become even smaller than the results with |AC| < 0.05. The 
situation is less serious for smaller value of #13, i.e., #13 ^ 3°, for which the correlation (5, 
C) is not so strong. It should be noted that we have assumed in our analysis that the 
detection efficiency does not decrease down to the neutrino energy E v ~ a few GeV, so if 
this assumption is not satisfied then the optimum muon energy may not be as low as Fig. 8 
indicates. 

In Figs. 7 and 8 we have taken 5 = tt/2 as a reference value. It is possible to do the 
same analysis for a value of 5 other than ir/2. The results for #13= 8°, 5° and 1° are given in 
Fig. 9 for three sets of the parameters (£ M =50GeV, L=3000km), (£ M =20GeV, L=1000km) 
and (£ M =20GeV, L=2000km). We observe that (^=50GeV, L=3000km) is better than 
(S M =20GeV, L=1000km) for smaller values of # 13 , but for larger values of 9 13 (E fl =20GeV, 
L= 1000km) can be more advantageous than the other. It should be emphasized that in all 
cases in Fig. 9 we can distinguish the case of 5 = tc from that of 5 = 0, since the necessary 
data size to reject S = is finite even for 5 = ir. This is because there are both contributions 
from sin 5 and cos 5 for the muon energy ^50 GeV. As we will see in section 4, for 
extremely high energy E^ 3>50 GeV we can show analytically that our A% 2 (CPV) becomes 
proportional to sin 2 S and distinction between 5 = n and 5 = is no longer possible. 

Sato et al. have been emphasizing the importance of the low energy option (E^ < 
lGeV), since the shape of sine curves in the CP violating probabilities becomes more con- 
spicuous for low energies. As can be seen from Fig. 7, the data set which is require to 
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reject a hypothesis 5 = is much larger than the value for the high energy case. This is 
because the small boost factor (P^/m^) 2 means that the number of events is statistically 
insufficient even though the CP conserving probabilities become relatively large in this case. 
Furthermore, in obtaining Fig. 7, we have assumed that the detection efficiency does not 
change for low energy. In reality, however, it becomes more and more difficult to distinguish 
wrong sign muons from right sign ones at lower energy so that the background fraction is 
expected to increase as — > small. Therefore, we conclude that the low energy option 
offers no advantage over the high energy region as long as we have the same detector size. 



Sato et al. [32] gave another criticism that neutrino factories at high energy are looking 
at only terms which are proportional to cos 5. However, as we will see in section 4, we 
can show analytically in the high energy limit that our Ax 2 (CPV), when minimized with 
respect to one of the oscillation parameters, is proportional to the square J 2 of the Jarlskog 
factor and has the same behavior oc 1/E^ as Ax 2 in the case of T violation. Therefore a 
naive argument that the square of difference of the quantity at 5 and the one at 5 = must 
depend only on cos S does not necessarily apply. 



III. ANALYSIS OF T VIOLATION 



It is known that measurements of T violation at neutrino factories are very difficult 
since electrons and positrons are difficult to distinguish at high energy. Here we discuss 
the possibility of measurements of T violation assuming perfect polarized muon beams for 
simplicity, i.e., P M cos# in eqs. (3) and (4) in (jTOJ is assumed to be ±1 so that either the flux 
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where x = E u /E^, so that we can compare the numbers of events (3/2) Njiy^ — > 
u e ; P^cosO = —1) from in decays of polarized /i~ and Nj(u e — > i/^; P^cosO = 0) from v e 
in decays of unpolarized fi + . Ax 2 in this case is thus given by 
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and the systematic errors are taken into account only for the /i + beam. In the definition ([5]) 
we do not optimize Ax 2 (TV) with respect to 9 k i, ^■'m 2 Ml C, since the numerator is defined 
only by quantities which are measured by experiments and a non zero value of Ax 2 (TV) 
immediately indicates 5^0. The date size to reject a hypothesis with 5 = at 3<r using 
A% 2 (TV) can be calculated as in the case of Ax 2 (CPV) and is shown if Fig. 10 for 5 = n/2 
and #i3 = 8° assuming the best fit values for (#12, Am^) and (623, Am^). The number of 
degrees of freedom in this case is one, i.e., the magnitude of the Jarlskog factor J = J sin 5, 
so that the data size is in general smaller than the case of CP violation for which we have 
six degrees of freedom. In practice, however, obtaining P M cos# ~ ±1 is technically difficult 
and the possibility to use T violation seems to be challenging. 



IV. LOW AND HIGH ENERGY BEHAVIORS OF A X 2 



In this subsection we will show analytically that the sensitivity to CP and T violation 
decreases as — > small (E^ <C lOGeV) or E^ — > large (E^ 3> 50GeV). Throughout this 
paper we assume sin 2 2# 13 > 1CT 3 (9 13 > 1°) so that we are always in the atmospheric regime 
in the language of P3fl , i.e., sin 2 2# 13 / sin 2 29 12 ^> (Am^/Am^) 2 . In this subsection we will 



ignore the effects of backgrounds and systematic errors for simplicity. 

To examine significance of CP/T violation analytically, we introduce the following sim- 
plified quantities: 



where 



A X 2 (CPV) = min A X 2 (6 M , Am 2 ,, S, C; 6 H , Am 2 ,, S = 0,C), 
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is defined as in (|Tj), 



(8) 



(9) 



are the number of events ((a, (3) = (e, /i) or (fi, e); in the case of (a, (3) = (/z, e) we 
assume perfect polarization as in the previous section so that the number of events is given 
by the same definition (|9|)) and we have ignored effects of the backgrounds and systematic 
errors and correlations of errors for simplicity in this section. Also we will assume that the 
cross section is proportional to the neutrino energy E v for any E u , i.e., a = a E u . Strictly 



speaking this assumption is not accurate, but it is known |38[ that < a < a E u is satisfied 
for low energy E v <C lGeV, so our approximation is sufficient to give an upper bound on 
the value of A% 2 for low energy. 
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Let us first look at the low energy limit (E v <C lOGeV). In this case matter effects are 
negligible and the probability can be replaced by that in vacuum. Thus we have 
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7rL 2 m 2 
(P(u e ^^;S)) 
12sl 3 sin 2 26 13 Ef t a 
irL 2 m 2 

h 1 



dx x (1 — x) sin 



V ^E u 



sm 



dx x (1 — x) sin 



2 ( Am 2 31 V 
\ ^E, , 



' Am 31 V 
AxE„. 



sm 



' Am 32 V 
AxE„ 



(12) 



(13) 



(14) 



where 2 = E v /E i 

sin 2 2^13/ sin 2 2^12 3> (Am^/Am^) 2 , and we have put #h 



we have assumed conditions for the atmospheric regime 

6 k f, Am 2 M = Amh, C 



C 



in fljjD instead of optimizing (P(f e — >■ ^v;<5)) — (-P(^e — * 8 — 0)) with respect to these 
variables, as that is sufficient to demonstrate that Ax 2 (CPV) decreases as E^ — > 0. If we 



keep L/E/j, fixed while L,E^ — > small, then all the quantities ([12]) , (O) and ( |I1 ) behave 
as so Ax 2 (CPV) oc £ M and A^ 2 (TV) cx E^ as ^ -»• with L/E^ fixed. Thus 

sensitivity to CP/T violation is asymptotically lost as E^ — > 0. This is consistent with our 
numerical results in previous sections. 

Next let us discuss the behavior of A% 2 in the high energy limit (E^ ^>50GeV). In this 
case we have to take into account the matter effect and we use the probability which has been 



obtained in J22[ to second order in O(0 13 ), 0(AE 21 /AE 31 ), 0(AE 21 /A) and 0(AE 21 L): 



P{v e -> v^S) ~ s^ 3 sin 2 2^ 13 



( AE- 



31 



V B 



sm 



>(f) 



2 2 n f AE 2 \ 

c 23 sin 26»i2 I 



sm 



/AL 
\~2~ 



+ 



~A£ 21 A£ 3 i . /AZA . / A£ 31 L\ 

8 J— — sm j sm j cos (* + j . 



(15) 



Since we assume sin 2 2^13/ sin 2 2#i2 ^> (Am^/ Am^) 2 = (AE21/AE31) 2 here, we can ignore 
the second term in fllBD . 
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It is straightforward to get the following high energy limit of Ax 2 (TV). Using ( |T5|) we 
have 



P(v e -> 5) - Piy^ -> u e ; 5) = P(v e -> S) - P(v e -> v^, -$) 

/ AT,\ 

(16) 



-Ami, (Ami) 2 L 9 / AL\ 

■ 2J e?, T* sm {!r) 



where we have expanded sm(AE 2 iL/2) ~ AE 2 \L/2 and have used the fact B = 
[(AE 31 cos 26»i3 - A) 2 + (AE 3 i sin26 l i3) 2 ] 1/2 - A as E u large. Therefore the number 
of events is given by 

(P{u e ^vs,S))-{P(v^v e ;S)) 

,2 /A_2 \2 



2AN a Jsm5Am 2 1 (Am 2 1 ) 1 2 / AL\ r 

vrm 2 A 2 L { 2 J J 1 %) 

** m , (17) 



12N a Jsm5Am 2 21 (Am 2 31 ) 2 1 . 2 (AL\ 
{P{u e ^v^8)) 



W y4 2 L 2 Sin2 I" J • (18) 



AW23 sin 2 2# 13 £^ . 2 /AL\ 

- - — Rin I I 

Hence we have the behaviors 

,2 \2 /a 2 \4 



A 2 mn iVoao288sin 2 ^J 2 (Am 2 1 ) 2 (Am 2 1 ) 4 1 . 2 /AL 
AX UVj " vrm 2 S 2 3 sin 2 2^ 13 E,A 2 ^ \ 2 

as — > large. 



(19) 



(|19|) indicates that the sensitivity to T violation decreases as E^ becomes very large. Also for 
a fixed large E^ A^ 2 (TV) is optimized for L ~ tc/A ~ 3x2000km/(p/2.7g-cm- 3 ) ~ 5000km. 
From numerical calculations we see that Ax 2 (TV) is optimized for (L, E^) ~ (3000km, 
50GeV) (see Fig. 1), so our analytic treatment is consistent with numerical calculations 
qualitatively. 

The behavior of Ax 2 (CPV) is a little more complicated, as we have to optimize A% 2 
with respect to 9 M , Am\ t , C. If we put Oki = Oki, Am 2 M = Am 2 ke , C = C as we did in (|12|). 
we have 

P{u e - v,- 6) - P{u e * Vfi'i 5 = 0) — 2J(coe* - l) Am ^f m 3i -L sin 2 (^) , (20) 

c- 4iV croJ(cos5 - llAm^Am?! K, . , /AL\ 
(P(i/ e - 6)) - (P(p. -,/„;* = 0)) = — —J- ^ 31 sm 2 J , 

ft 

(21) 

so that we naively have the following behavior 
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Ax 2 (naive CPV) 



N a 32(cos<5 - if. J 2 (Am 2 21 ) 2 [Am^f E, 



Tim, 



s 2 23 sin 2 20 13 



31 ) ^fi 

A 2 L 2 8 V 2 



2 (AL 



(22) 



It turns out that it is sufficient to consider the correlation of two variables (5, X), where 
X is 9ke, Am 2 ^ or C, to demonstrate Ax 2 (CPV) oc l/E^. Except for the correlations (5, 
0i 2 ) and (5, Amfi), we can ignore terms of order 0({AE2\/ AE 3 i) 2 ). From the assumption 
sin 2 20 13 / sin 2 20 12 3> (Am^/Am^) 2 , (|T5| ) is approximately given by 



. AL31 . ML 

S93 sm 20 13 — - — sin 

3 A V 2 



+ 



4J A# 



21 



(") 



sm I — — I cos 5 



AE31L 



(23) 



s 23 sin 2^13 A 

where we have used A — AE 3 i ~ A for E v — ► large, and we have ignored terms of order 
C((Ai?2i/A£^3i) 2 ). In the case of the two variable correlation (5, #13), to minimize the 
square of 

P(u e -> z/ M ; 6»i3, <5) - P(z/ e -> 9 13 , 5) 
AE 31 . [AL\ A J 



g 2 3 sin 2fli3 ^ sin ( — j 

* AE31 
s 23 sm20i 3 



A£ 21 . /AL\ / AE 31 L 

, , — - — sm cos H 

V 2 7 s 23 sin20 13 A V 2 / V 2 



sm {-TJ + 



4 J A£ 21 . (AL\ AE 31 L 

: — sm cos H 

s 23 sin20 13 A \ 2 J \ 2 



it is sufficient to takeQ 



4J 



Am 2 x 



sin 2^i3 = sin 26 13 - 

s 23 sm 2^13 Am 31 



(cos 5 — cos 5 



(24) 



(25) 



where we have used in (p4j) and ( |25|) the fact J/ sin 20x3 = cos0i3 x (independent of #13) ~ 
cos 0i3 x (independent of i3 ) which holds because sin 2 13 <C 1. Notice that the phase 
AE 31 L/2 which appears together with 5 in cosine in (|24j ) disappears as Lj, — > large. Plugging 
(BSD in (Mi, we find 



2 



V£ 013, 5) - P&e 

2 sin 2 (^) (sin 2 20 13 



31 



+ 8J 



23 V a 

~AE 21 AE 31 . 2 /AL 
/ sm I 



A 2 



V 2 



cos f 5 + 



- sin 2 20 
AE 31 L 



13 



COS I 5 



AE31L 



~ AE21 AE 3 \ 
8J A 2 Sm V 2 



2 /AL 



cos ( 5 



AE 31 L 



U 



~AE2\AE 3 i 



A 2 



sm 



— cos ( 5 
(AL\ 



AE 31 L 



sin 5 — sin 5 



— cos 5 + cos 5 
AE 3i L 



(26) 



3 Here we do not discuss the other solution of the quadratic equation which was discussed by [44], 
since we are mainly interested in rejecting 5 = rather than determining the precise value of 5. 
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where we have expanded sin(A£ 3 iL/2) ~ AE 31 L/2, cos(AE 31 L/2) - 1 ~ - (AE 31 L) 2 /2 ~ 
in the last step in (p6|) . Hence we get 

{P(u e -> ^; 13 , 5) - P(z/ e -> 13 , 5 = 0)) 
= ^^K-tM. (27, 



We see from (|27| ) that if we optimize Ax 2 (CPV) with respect only to #i 3 then Ax 2 (CPV) 
behaves as 

vrm 2 s^ 3 sm 26»i 3 E^A 2 

Note that the behavior of Ax 2 (CPV; (5, #13)) which is optimized with respect to #i 3 is 
quite different from that of A% 2 (naive CPV) in fl2"2"|). We observe that the dependence 
of Ax 2 (CPV; (5, #13)) on is the same as that of Ax 2 (TV). It should be also empha- 
sized that Ax 2 (CPV; (6,613)) is proportional to sin 2 6 and does not depend on cos 6 unlike 
Ax 2 (naive CPV) in ©. 



We can play the same game for 623, Am§ 2 and C. In the case of the two variable 
correlation (6, 623), 

A: J A/72 2 

sin 623 = sin 623 - 2 . f (cos 6 - cos 6) (29) 

s 2 3 sin 2tfi 3 Am 31 v / 

minimizes Ax 2 (CPV; (6,62s)) an d we have 

A^(CPV; (6, ft,)) g ^ S1 " 2 ^ (AmIl) ' Sm ;'- 4 f/ 2) ■ (30) 

vrm 2 s^sm 20i3 EpA 2 

which is the same as Ax 2 (CPV; (6, 613)). In the case of the two variable correlation (6, 
Am| 2 ), using 



P(v e -> ^; 5) 

we find 



A£ 31 (AL\ 4 J A£ 21 ( AL- 1 

-I sin 



S23 sin 26*13 — ; — sin H — sin I cos 6 

A V 2 J s 23 sin20i 3 A V 2 1 



(31) 



4j , _ . 

Am 3 i = Am^ 2 — ~ — ^— Am 21 (cos 5 — cos 5) (32) 

s 23 sin 26*13 

minimizes Ax 2 (CPV; (6, Am 31 )). We obtain 

A/(CPV; (d, H„ = ^ ^^^(A^)^W) 3 

Trm 2 s| 3 sm 20 i3 E^A 2 
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which again is the same as Ax 2 (CPV; (5, #13)). In the case of the two variable correlation 



sin(AL/2) _ s in(AL/2) 



A 



A 



1 - 



4J Amli / ? 
~~9 — . 9 ttt, — I — 5- cos d — cos 
s 2 23 sin' 20 i3 Amfi V 



(34) 



minimizes Ax 2 (CPV; (<5, C)) and we get 



18N a J 2 sin 2 5 (Am^) 2 (Am^ sin 2 (AL/2) 



,2 \2, 



A X 2 (CPV; (5,C)) 



7rmf 



s| 3 sin 2 20i3 



i^ 2 



which once again is the same as Ax 2 (CPV; (<5, (9 13 )). The expressions (E3| 



(35) 



31) and 



([34]) for the optimal values for #13, #23, Am^ and C explain why the correlation has a cosine 
curve for large and small L in Figs. 1, 2, 4 and 6. 

In the case of the correlations (5, 612), and (5, Am^), we have to take into account of 
terms of order 0((AE 2 i/ AE 3 i) 2 ) . For (5, 0\ 2 ), we have 



. - 4 J Am 2 - 

sin 20i2 = — 7T- -7 — 5- cos + 
c{ 3 A< 



f«^f) cos , J + WAn4 
\ V c 23 Am^y q 3 Am^ 



cos 5 sin 26 12 + sin 2 2(9 12 , (36) 



and this optimizes Ax 2 (CPV). We find 



012, <5)-P(l/e-> ^012, 



8J A£ 2 iA£ 3 i . 2 /AL\ 



sin 20 
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A 2 Sm V 2 / 



sin 20x2 cos ( 5 + 



AE 31 L 



AE 31 L 

sin 26*12 cos ( 5 H — ) — sin 20 i2 cos 5 + sin 29 12 cos 5 



: — sm I sm sm 20i 2 — sm sm 20 u 

in 9flio A 2 \ 2 / ^ 



AE 3 iL 



sm20 12 A 2 V 2 / V 2 

where we have expanded sin(Ai?3iL/2) ~ AE 3 \Lj2. By putting 5 = 0, we obtain 

18iYo<x J 2 sin 2 5 (Am 2 x ) 2 (Am^) 2 sin 2 {AL/2) 



(37) 



A X 2 (CPV; (c5,(9 12 )) 



7rmf 



s| 3 sin 2 20 i3 



£^ 2 



(38) 



which once again is the same as Ax (CPV; (5, 6*13)). For (5, Am 21 ), we have 



4J 



c{ 3 sin' 20 



Am^ cos 5 



12 



- \ I 2 ^ofl ) cos 2 ^(Am| 1 ) 2 + 2 8 { Q . Am 2 1 Am 2 1 cos(5 + (Am 2 1 ) 2 (39) 
\ V c 23sm 2012/ c^sm 20 i2 



which leads to 



P{y e -> Am 21 , 5) - P(^ e -> z^i Am\ x , 5) 



p ~ AE 21 AE 31 . 2 {AL\ ( - 

8J --Tj sm I — I (sino -smo 



AE 31 L 



(40) 
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Thus we get 



18N a J 2 sin 2 S (Am 2 2l f (Am 2 31 ) 2 sin 2 (AL/2) 



A X 2 (CPV; (5, Amy) * — ^ f' / ; . (41) 

which once again is the same as Ax 2 (CPV; (<5, #i 3 )). Unlike the cases for (5, #i 3 ), (5, 623), 
(<5, Am^) and (5, C), the optimal values (|36D and (|39| ) have nontrivial behaviors even for 
large and small L, as we can see from Figs. 3 and 5. 

We have seen analytically that two variable correlations give us the behavior 
Ax 2 (CPV) oc sm 2 5/E^ and this behavior is the same as Ax 2 (TV). Although it is diffi- 
cult to discuss correlations of more than two variables analytically, the discussions above are 
sufficient to demonstrate that sensitivity to CP violation decreases as becomes larger. 
In fact we have verified numerically that Ax 2 (CPV) decreases as the muon energy increases 
(E^ > lOOGeV). The conclusion in this subsection is qualitatively consistent with the work 
|36| by Lipari who claims that sensitivity to CP violation decreases as becomes large. 
However it may not be quantitatively consistent with |JB| in which it was suggests that 



sensitivity starts getting lost for E u > a few GeV. In our discussion here it was necessary 
to have |Ai?3iL| <C 1 which may not be attained for L ~ 3000km and E^ < 50GeV. Our 
numerical calculations in the previous section indicate that the sensitivity is optimized for 
20GeV ^ < 50GeV which is quantitatively consistent with the results in [|2T| |3TH . This 
interval for is the intermediate energy region which cannot be treated analytically using 
our arguments in this section. In fact it seems difficult to explain analytically the strong 
correlation of (6, C) for E^ ~50 GeV and L ~ 3000km. (cf. Fig. 1) 



V. JHF EXPERIMENT 



The JHF project [18| has been proposed to perform precise measurements of the oscil- 



lation parameters. The possible extension of this project includes the upgrade of the power 
to 4MW and the construction of a mega-ton detector |[39|| . The possibility to measure CP 



violation at the JHF project has been discussed by [|45| , [l2| , [46[ . Here we briefly discuss the 
possibility of measurements of CP violation at the JHF experiment with power 4MW and 
a 1 Mton detector as a comparison with neutrino factories. As in previous sections, we will 
take into consideration the correlations of all the oscillation parameters. In the case of the 
JHF experiment, which has the baseline L ~ 300km, the matter effect is almost negligible 
and it is possible to compare the numbers of events for — > v e and — > v e directly by tak- 
ing into account the difference of the cross sections between <j vN and o vN . However, we use 
the same A% 2 as in section 2, as discussions with the same criterion gives more transparent 
comparisons between neutrino factories and the superbeam at JHF. 

The correlations of two variables (5, X), where X is 8ke, Amy or C, are shown in 
Fig. 11, where the central values for these parameters are those of the best fit point, i.e., 
sin 2 2^i 2 = 0.75, Am 2 21 = 3.2 x 10" 5 eV 2 ; sin 2 2 6 23 = 1.0, Am 2 2 = 3.2 x 10~ 3 eV 2 , C=1.0 and 
we have used a reference value #13 = 8°. In this calculation the narrow band beam (NBB) 
(the flux referred to as LE27T in [|18|]) is used, and it is assumed for simplicity that there are 
no backgrounds and the detection efficiency is 70% in Fig. 11. Note that for the purpose of 
measurements of CP violation NBB is more advantageous than the wide band beam, as the 
former has better energy resolution. 
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As in the section 2, we have evaluated numerically the data size required to reject a 
hypothesis with 5 = 0. Of course the data size depends on the true value 5 and the results 
obtained by varying the six variables (9^, Am|£, C) are plotted in Fig. 12, where we have 
taken the best fit values for (0 12 , Am 2 21 ), (0 23 , Am|J, 13 = 8°, 5°, 1°(2°), and the NBB is 
used. The vertical axis of Fig. 12 stands for the data size required per ktx(z/ M 1 year + 
Up 2 years). We have used two ways of v e selections, one is 1-ring e-like selection which 
has the background fraction /# = 1.8 x 10~ 2 , the detection efficiency 70.4%, and the other 
one is n° cut selection which has the background fraction /# = 2 x 10~ 3 , the detection 
efficiency 50.4% ]4"5| . In the case of the 1-ring e-like selection, for 6* 13 = 1° the systematic 
error becomes so large that the data size required to reject 5 = becomes infinite. Also in 
this case the number of events for 5 = n becomes almost the same as that for 5 = up to 
the systematic errors and there is no way to distinguish the case of 5 = ir and that of 5 = 0. 
However, as long as the value of 5 is not close to or ir and 13 > 3°, the JHF with 4MW 
power and a 1 mega ton detector will be able to demonstrate 5 ^ at 3crCL. 



VI. DISCUSSIONS 



In order to be more concrete, we need the knowledge on the uncertainty of the matter 
effect A. The error of A = \^2GpY e p comes from those of Y e and p. The error of Y e has 
been discussed by j|7| and it is about 2% and geophysicists |48|j49| agree with it. Without 



any uncertainty of the matter effect, it has been claimed that a medium baseline experiment 
(L ~ 3000km, ~ 50GeV) is best for measurements of CP violation. In that case the depth 
of the neutrino path is at most 200km and most of the neutrino path is in the upper mantle. 
It is known in geophysics |5(J that the crust has relatively large latitude-longitude dependent 



fluctuations around constant density. On the other hand, in the case of the upper mantle, 



some geophysicists claim that fluctuations around constant density are a few % [51,48,49 



while another |52|] says that they may be as large as 5 %. However, such discussions are 



based on normal mode studies in seismology which are confined to long wavelength features, 



and it was pointed out [Hj,^ that the fluctuations in the density in the analysis of neutrino 
factories may be larger than 5 %, since the width of the neutrino beam is much smaller than 
typical wavelengths in seismological studies. If that is the case, then it follows from Fig. 8 
that the case L ~1000km is better than the case L ~3000km, since the former is insensitive 
to the uncertainty of the matter effect. More detailed analysis with seismological discussions 
seems to be necessary to determine the optimal baseline and the muon energy in neutrino 
factories. 

We have also seen that if there is a nontrivial fraction of backgrounds, which is expected 
in the case of water Cherenkov or liquid argon detectors, then medium energy (lOGeV 
;$ EptSi 20GeV) is more advantageous than high energy (E^ ~ 50GeV), as systematic errors 
become so large in the latter case. 

The bottom line of the present paper is that either the high or medium energy option 
is certainly better than the low energy option which has been advocated by some people 
|32|,|34|]. We have arrived at this conclusion on the assumption that the energy threshold 
is as low as 0.1 GeV, and the detection efficiency is independent of the neutrino energy. 
In practice, it may be very difficult to have such a low threshold and to keep such a good 
detection efficiency down to 0.1 GeV, so it is expected that the low energy option (E, 



< 
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lOGeV) becomes less and less advantageous. 

If we can realize a nearly perfectly polarized muon beam, it may be possible to measure 
T violation, which is much better than CP violation, since we do not have to worry about 
the uncertainty of the density of the Earth or other oscillation parameters and a nonzero 
value of T violation immediately indicates 5^0. Efforts should be made along these lines. 

If #13 ~ 3° and if the value of 5 is not close to or it, then the JHF experiment with 
4MW power and a 1 mega ton detector will be able to demonstrate 5 ^ at 3cxCL. On 
the other hand, if #13 < 3°, then neutrino factories seem to be the only experiment which 
can demonstrate 5^0. In that case, depending on the situation such as the fraction of 
backgrounds, the uncertainty of the matter effect and the magnitude of #13, the option with 
(Efj, ~50GeV, L ~3000km) may be advantageous (or disadvantageous) over (E^ ~20GeV, 
L ~1000km). In both a neutrino factory with (20GeV< <50GeV, 1000km< L <3000km) 
and the JHF experiment, our A% 2 (CPV) depends not only on sin 5 but also cos 5, [] so that 
we can in principle distinguish 5 = tt from 8 = as long as the statistical significance 
overcomes the systematic errors. This is not the case for a neutrino factory with large 
systematic errors for small #13, i.e., for i? M =50GeV and #13 = 1° (cf. Fig. 9), and for the 
JHF experiment with less S/N ratio i.e., when the 1-ring e-like selection is adopted, or when 
#13 = 8° and the 7r° cut selection is adopted (cf. Fig. 12). 

In this paper we adopted simplified assumptions such as that the detection efficiency is 
independent of the neutrino energy, that the threshold energy can be taken as low as O.lGeV, 
and that the uncertainty of the matter effect is at most 5 %, but we need much more detailed 
experimental information to obtain the optimum muon energy and the baseline. 



NOTE ADDED 



Just before we hit a return key to submit our paper to the preprint archive, we became 
aware of the work by M. Freund et al. 



53fl , where the similar topics has been discussed from 



a slightly different viewpoint. 
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4 If we evaluate Ax 2 (CPV; (5, #13)) in (p8|) to the next leading order in Am^L/E^ then J 2 sin 2 S in 
(pS[) is replaced by J 2 (sin5+const.(Am 2 1 L/£' At ) cos <5) 2 , and Am 2 1 L/£' M is not necessarily negligible 
in either case. 
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Figures 

Fig. la Correlations of errors of 5 and the normalization C for L=100km, 1000km, 2500km, 
6300km and for £ M =3GeV, 20GeV, 50GeV. A X 2 = 0.18,0.37,0.73 corresponds to 
lcCL, 90%CL, 99%CL, respectively for two degrees of freedom. The oscillation 
parameters are Am^ = 1.8 x 10~ 5 eV 2 , Am^ = 3.5 x 10~ 5 eV 2 , sin 2 29n = 0.76, 
sin 2 2^23 = 1.0, #i3 = 8°, 5 = ir/2. The number of useful muon decays is 10 21 /i ■ lOkt. 
No backgrounds are taken into consideration in Figs. 1-6. 

Fig. lb The same correlation as Fig. la for #13 = 5°, 1°. The oscillation parameters and 
other reference values are the same as in Fig. la. 



Fig. 2 Correlations of errors of 5 and #13. 
values are the same as in Fig. la. 

Fig. 3 Correlations of errors of 5 and de- 
values are the same as in Fig. la. 

Fig. 4 Correlations of errors of 5 and 623- 
values are the same as in Fig. la. 



The oscillation parameters and other reference 
The oscillation parameters and other reference 
The oscillation parameters and other reference 



Fig. 5 Correlations of errors of S and Am^. The oscillation parameters and other reference 
values are the same as in Fig. la. 

Fig. 6 Correlations of errors of 5 and Am|. The oscillation parameters and other reference 
values are the same as in Fig. la. 

Fig. 7 The contour plot of equi-number of data size required (in the unit of kt) to reject a 
hypothesis 5 = at 3a using Ax 2 (CPV) (|3|) in the case of a neutrino factory with 10 21 
useful muon decays, the background fraction f B = 10~ 5 or 10~ 3 , 9 13 = 8°, 5°, 1°. The 
other oscillation parameters are the same as in Fig. 1. 

Fig. 8 The same as Fig. 7 with #13 = 8°, except that the uncertainty of the matter effect is 
assumed to be larger \AC\ < 0.1 or |AC| < 0.2. 

Fig. 9 The number of data size required to reject a hypothesis 5 = at 3cr for a neutrino 
factory using Ax 2 (CPV) @ as a function of the true value of 5 for fg = 10~ 3 . All 
the assumptions except for S are the same as in Fig. 7. The situation is improved for 
smaller # 13 . 

Fig. 10 The contour plot of equi-number of data size required (in the unit of kt) to reject 
a hypothesis 5 = at 3cr using T violation (J^) in the case of a neutrino factory with 
10 21 useful muon decays, the background fraction fg = 10~ 5 or 10~ 3 , #13 = 8°, 5°, 
1°. The other oscillation parameters are the same as in Fig. 1. Crucial assumption in 
this case is that the polarization is perfect, and under this assumption T violation is 
3 times better than the CP violation analysis in section 2. 
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Fig. 11 The correlations of errors of (5, C), (5, 9 ri ), (5, #i 2 ), (5, 623), (<5, Am^), (5, Ami) in 
the case of the JHF experiment with 4MW power, a 1 mega ton detector and NBB. No 
backgrounds are taken into consideration in these figures. The oscillation parameters 
used are the same as in Fig. 1. 

Fig. 12 The number of data size required (in the unit of kt) to reject a hypothesis 5 = at 
3cr for the JHF experiment with 4MW power, a 1 mega ton detector and NBB using 
A% 2 (CPV) (|3|) as a function of the true value of 5. Unlike in the case of Fig. 11, the 
effects of backgrounds are taken into account in this figure. The oscillation parameters 
used are the same as in Fig. 1, and two ways of cuts (1-ring e-like and tt° cut) |55| 



are used. In the case of the 1-ring e-like selection, 6*13 = 1° does not have a solution 
because the systematic errors become so large. 
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